Strategies for wheat stripe rust pathogenicity identified by transcriptome sequencing

Stripe rust caused by the fungus Puccinia striiformis f.sp. tritici (Pst) is a major constraint to wheat production worldwide. The molecular events that underlie Pst pathogenicity are largely unknown. Like all rusts, Pst creates a specialized cellular structure within host cells called the haustorium to obtain nutrients from wheat, and to secrete pathogenicity factors called effector proteins. We purified Pst haustoria and used next-generation sequencing platforms to assemble the haustorial transcriptome as well as the transcriptome of germinated spores. 12,282 transcripts were assembled from 454-pyrosequencing data and used as reference for digital gene expression analysis to compare the germinated uredinospores and haustoria transcriptomes based on Illumina RNAseq data. More than 400 genes encoding secreted proteins which constitute candidate effectors were identified from the haustorial transcriptome, with two thirds of these up-regulated in this tissue compared to germinated spores. RT-PCR analysis confirmed the expression patterns of 94 effector candidates. The analysis also revealed that spores rely mainly on stored energy reserves for growth and development, while haustoria take up host nutrients for massive energy production for biosynthetic pathways and the ultimate production of spores. Together, these studies substantially increase our knowledge of potential Pst effectors and provide new insights into the pathogenic strategies of this important organism.J.P.R. is an ARC Future Fellow (FT0992129). This project has been supported by Bioplatforms Australia through funding from the Commonwealth Government NCRIS and Education Investment Fund Super Science programs

Crystallization, X-ray diffraction analysis and preliminary structure determination of the TIR domain from the flax resistance protein L6

The Toll/interleukin-1 receptor (TIR) domain is a protein-protein interaction domain that is found in both animal and plant immune receptors. In animal Toll-like receptor signalling, both homotypic TIR-domain interactions between two receptor molecules and heterotypic interactions between receptors and TIR-domain-containing adaptors are required for initiation of an innate immune response. The TIR domains in cytoplasmic nucleotide-binding/leucine-rich repeat (NB-LRR) plant disease-resistance proteins are not as well characterized, but recent studies have suggested a role in defence signalling. In this study, the crystallization, X-ray diffraction analysis and preliminary structure determination of the TIR domain from the flax resistance protein L6 (L6TIR) are reported. Plate-like crystals of L6TIR were obtained using PEG 200 as a precipitant and diffracted X-rays to 2.3 angstrom resolution. Pseudo-translation complicated the initial assignment of the crystal symmetry, which was ultimately found to correspond to space group P2(1)2(1)2 with two molecules per asymmetric unit. The structure of L6TIR was solved by molecular replacement using the structure of the TIR-domain-containing protein AT1G72930 from Arabidopsis as a template

The ins and outs of rust haustoria

Rust diseases caused by fungi of the order Pucciniales afflict a wide range of plants, including cereals, legumes, ornamentals, and fruit trees, and pose a serious threat to cropping systems and global food security. The obligate parasitic lifestyle of these fungi and their complex life cycles, often involving alternate hosts for the sexual and asexual stages, also make this group of pathogens of great biological interest. One of the most remarkable adaptations of rust fungi is the specialized infection structure that underpins the sustained biotrophic association with hosts; the haustorium.This work has been supported by grants from the TwoBlades corporation and the Australian research council

The genome sequence and effector complement of the flax rust pathogen Melampsora lini

Rust fungi cause serious yield reductions on crops, including wheat, barley, soybean, coffee, and represent real threats to global food security. Of these fungi, the flax rust pathogen Melampsora lini has been developed most extensively over the past 80 years as a model to understand the molecular mechanisms that underpin pathogenesis. During infection, M. lini secretes virulence effectors to promote disease. The number of these effectors, their function and their degree of conservation across rust fungal species is unknown. To assess this, we sequenced and assembled de novo the genome of M. lini isolate CH5 into 21,130 scaffolds spanning 189 Mbp (scaffold N50 of 31 kbp). Global analysis of the DNA sequence revealed that repetitive elements, primarily retrotransposons, make up at least 45% of the genome. Using ab initio predictions, transcriptome data and homology searches, we identified 16,271 putative protein-coding genes. An analysis pipeline was then implemented to predict the effector complement of M. lini and compare it to that of the poplar rust, wheat stem rust and wheat stripe rust pathogens to identify conserved and species-specific effector candidates. Previous knowledge of four cloned M. lini avirulence effector proteins and two basidiomycete effectors was used to optimize parameters of the effector prediction pipeline. Markov clustering based on sequence similarity was performed to group effector candidates from all four rust pathogens. Clusters containing at least one member from M. lini were further analyzed and prioritized based on features including expression in isolated haustoria and infected leaf tissue and conservation across rust species. Herein, we describe 200 of 940 clusters that ranked highest on our priority list, representing 725 flax rust candidate effectors. Our findings on this important model rust species provide insight into how effectors of rust fungi are conserved across species and how they may act to promote infection on their hosts.This work was funded by a grant from the CSIRO Transformational Biology Capability Platform to Adnane Nemri. Claire Anderson was supported by an ARC Discovery Grant (DP120104044) awarded to David A. Jones and Peter N. Dodds

Crystallization and X-ray diffraction analysis of the C-terminal domain of the flax rust effector protein AvrM

The flax rust effector AvrM is a secreted protein of unknown fold that is recognized by the M resistance protein in flax. In order to investigate the structural basis of the AvrMM interaction and possible virulence-associated functions of AvrM, the C-terminal domains of two different AvrM variants (AvrM-A and avrM) were crystallized. Crystals of native AvrM-A were obtained using pentaerythritol ethoxylate (15/4 EO/OH) as a precipitant and diffracted X-rays to 2.9 angstrom resolution. Selenomethionine-derivative crystals of similar quality were obtained using PEG 1500 as a precipitant. Both the native and selenomethionine-labelled AvrM-A crystals had symmetry of space group C2221 with eight molecules in the asymmetric unit. Crystals of avrM had symmetry of space group P212121 and diffracted X-rays to 2.7 angstrom resolution. Initial AvrM-A phases were calculated using the single-wavelength anomalous dispersion (SAD) method and a partial model was built. Phases for avrM were obtained by molecular replacement using the partial AvrM-A model

Comparative genomics of Australian isolates of the wheat stem rust pathogen Puccinia graminis f. sp. tritici reveals extensive polymorphism in candidate effector genes

The wheat stem rust fungus Puccinia graminis f. sp. tritici (Pgt) is one of the most destructive pathogens of wheat. In this study, a draft genome was built for a founder Australian Pgt isolate of pathotype (pt.) 21-0 (collected in 1954) by next generation DNA sequencing. A combination of reference-based assembly using the genome of the previously sequenced American Pgt isolate CDL 75-36-700-3 (p7a) and de novo assembly were performed resulting in a 92 Mbp reference genome for Pgt isolate 21-0. Approximately 13 Mbp of de novo assembled sequence in this genome is not present in the p7a reference assembly. This novel sequence is not specific to 21-0 as it is also present in three other Pgt rust isolates of independent origin. The new reference genome was subsequently used to build a pan-genome based on five Australian Pgt isolates. Transcriptomes from germinated urediniospores and haustoria were separately assembled for pt. 21-0 and comparison of gene expression profiles showed differential expression in ∼10% of the genes each in germinated spores and haustoria. A total of 1,924 secreted proteins were predicted from the 21-0 transcriptome, of which 520 were classified as haustorial secreted proteins (HSPs). Comparison of 21-0 with two presumed clonal field derivatives of this lineage (collected in 1982 and 1984) that had evolved virulence on four additional resistance genes (Sr5, Sr11, Sr27, SrSatu) identified mutations in 25 HSP effector candidates. Some of these mutations could explain their novel virulence phenotypes.Authors wish to thank the Two Blades Foundation for financial support. Part of this work was supported through access to facilities managed by Bioplatforms Australia and funded by the Australian Government National Collaborative Research Infrastructure Strategy and Education Investment Fund Super Science Initiative

Geometry of River Networks II: Distributions of Component Size and Number

The structure of a river network may be seen as a discrete set of nested sub-networks built out of individual stream segments. These network components are assigned an integral stream order via a hierarchical and discrete ordering method. Exponential relationships, known as Horton's laws, between stream order and ensemble-averaged quantities pertaining to network components are observed. We extend these observations to incorporate fluctuations and all higher moments by developing functional relationships between distributions. The relationships determined are drawn from a combination of theoretical analysis, analysis of real river networks including the Mississippi, Amazon and Nile, and numerical simulations on a model of directed, random networks. Underlying distributions of stream segment lengths are identified as exponential. Combinations of these distributions form single-humped distributions with exponential tails, the sums of which are in turn shown to give power law distributions of stream lengths. Distributions of basin area and stream segment frequency are also addressed. The calculations identify a single length-scale as a measure of size fluctuations in network components. This article is the second in a series of three addressing the geometry of river networks.Comment: 16 pages, 13 figures, 4 tables, Revtex4, submitted to PR